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Quantum phonon transport of molecular junctions amide- linked with carbon 

nanotubes: a first-principle study 
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Quantum phonon transport through benzene and alkane chains amide-linked with single wall car- 
bon nanotubes (SWCNTs) is studied within the level of density functional theory. The force constant 
matrices are obtained from standard quantum chemistry software. The phonon transmission and 
thermal conductance are from the nonequilibrium Green's function and the mode- matching method. 
We find that the ballistic thermal conductance is not sensitive to the compression or stretching of 
the molecular junction. The terminating groups of the SWCNTs at the cutting edges only influence 
the thermal conductance quantitatively. The conductance of the benzene and alkane chains shows 
large difference. Analysis of the transmission spectrum shows that (i) the low temperature thermal 
conductance is mainly contributed by the SWCNT transverse acoustic modes, (ii) the degenerate 
phonon modes show different transmission probability due to the presence of molecular junction, 
(iii) the SWCNT twisting mode can hardly be transmitted by the alkane chain. As a result, the 
ballistic thermal conductance of alkane chains is larger than that of benzene chains below 38 K, 
while it is smaller at higher temperature. 

PACS numbers: 05.60.Gg, 44.10.4-i, 63.22.-m, 65.80.-|-n 



I. INTRODUCTION 

The electronic transport properties of a molecular de- 
vice depend much on the underlying nuclear configura- 
tion and the electron-phonon couplingi. The mutual 
interaction between them induces Joule heating in the 
molecule, which may have important consequences on the 
functionality and reliability of a molecular device. Study 
of the heat transport by phonons in molecular junctions 
is crucial for a better understanding of their electronic 
transport properties. Phonon thermal transport itself is 
also interesting. Detailed understanding of the underly- 
ing transport mechanism is especially useful for the de- 
sign of novel phononic devices^. Furthermore, combined 
study of the electronic and phononic transport in molecu- 
lar junctions is the first step toward the design of molecu- 
lar thermoelectric devices^. Great experimental progress 
has been made in these directions recently, which enables 
one to measure the thermal and thermoelectric trans- 
port properties of molecular junctions^i^. Advances in 
experimental technique call for new theoretical method 
to predict the thermal conductance of molecular junc- 
tions. Although scmicmpirical or minimal model calcu- 
lation is helpful to understand the underlying physics^ii. 
For a detailed quantitative study, a parameter-free, first- 
principle method is highly desirable^. Furthermore, for 
many molecular structures there exists no empirical inter- 
atom potential. 

In this paper wc introduce such a method based 
on a standard quantum chemistry software, the Gaus- 
sianOS parckage^. Given the molecule structure, we 
can obtain the force constant matrices after perform- 
ing the energy minimization. The thermal conduc- 
tance can be calculated via available theoretical methods. 
Among them are the nonequilibrium Green's function 



(NEGF) mcthodiSiiiii^ii^ii^, which has been successfully 
used to predict the electronic conductance of molecular 
junctionsi^ii^iii. To study the transport of each phonon 
mode, we will also use the mode-matching methodic, 
which is equivalent to the NEGF method in the ballistic 
limii^i^. Using these methods, we first compare the ther- 
mal conductance of a benzene ring amide-linked with two 
(6, 0) single-walled carbon nanotubes (SWCNTs) under 
different compression or stretching conditions. Then we 
study the effect of SWCNT terminating group at the cut- 
ting edges. Finally, we compare the phonon transmission 
probability and the thermal conductance of benzene and 
alkane chains with the same leads. Although the electron 
contribution to the thermal conductance may be compa- 
rable with or larger than that of phonons, inclusion of 
this effect is out of the scope of present study. 



II. MOLECULAR STRUCTURE AND 
THEORETICAL METHOD 

In this section we first introduce the system Hamilto- 
nian and the procedure to obtain it from the GaussianOS 
package. Then we briefly outline the NEGF, the mode- 
matching method, and their relationship. For a detailed 
discussion, we refer the reader to Ref. [ij. The system we 
arc interested in is a molecular junction connected with 
two periodic semi-infinite leads at both sides. It is a stan- 
dard treatment to divide the whole structure into three 
parts: the center junction and two leads acting as ther- 
mal baths. The boundaries between the center and the 
baths may be at arbitrary positions, and not correspond 
to any physical interface. But it is desirable to include 
part of the leads into the center region, since we need to 
make sure that there is no direct interaction between the 



two baths, which is required by the NEGF formaUsm. 
By doing this, we can also include the charge transfer 
effect between the leads and the molecule junction. In 
this setup, the system Hamiltonian can be written as 



the formula) 



a=L,C,R 



^\'rv^^,f + {vPfv^^u'' + K, (1) 



where H^ = ^l""") ''^^ + ^l""") K°'u°' represents har- 
monic oscillators, u" is a column vector consisting of 
all the mass normalised displacement variables in region 
a, and u" is the corresponding conjugate momentum. 
K" is the spring constant matrix in the tight-binding 
form and V^'~^ — {V'~^^)^ is the coupling matrix of the 
left lead to the central region; similarly for V'-^^ . Vn 
is the nonlinear interaction in the center, which could 
be 1^ = ^ J2ijk tijkUiUjUk for cubic nonlinearality. We 
ignore the nonlinear interaction in this paper and only 
briefly discuss its effect in Sec. IIIII 

We study the phonon thermal conductance of ben- 
zene and alkane chains covalently bonded with two (6, 0) 
SWCNTs via the amide group. This is relevant to a re- 
cent experimental setup22i, where the SWCNT is oxida- 
tively cutted, and the cutting edges are covalently bonded 
with molecular chains via the amide group. To get the 
force constants of the system, we do two separate runs for 
the center and the leads using GaussianOS. For the cen- 
ter, we include extra one and a half periods of SWCNT at 
each side, which is terminated by hydrogen atoms at the 
outer boundaries. The cutting edges may be terminated 
by H or COOH group. We optimize the center at the 
level of b31yp density functional method using the 6-31G 
basis set. We first relax the structure by constraining 
the leads to be coaxial to get an optimized distance be- 
tween the two leads. Then we fix the outer half period of 
the SWCNTs and the H atoms at both sides and allow 
all other atoms to move freely. The optimized struc- 
ture is used to get the dynamical matrix. The inset of 
Fig. [2] shows the optimized benzene structure terminated 
by the H and COOH groups at the cutting edges. For 
the lead, it is desirable to use periodic boundary condi- 
tion to compute the dynamical matrix, while this is out 
of the ability of GaussianOS. So we optimise 7 SWCNT 
periods and extract the force constant from the central 
period to minimize the finite size effect. To connect the 
molecule with the two leads, we remove the outer fixed 
atoms in the molecule and connect the remaining period 
of SWCNT with one semi-infinite SWCNT at each side. 
For the coupling matrix between the center and the leads, 
we only include coupling between one period of SWCNT 
atoms in the center and one period in the leads. This 
is a good approximation so long as we include enough 
SWCNT atoms into the center molecule. 



In the NEGF method as described in Refs.|10||ll|UJll3|, 
Il4l . thermal conductance of the molecular junction can be 
calculated from the Landauer formula {h = I throughout 
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with the transmission coefficient 



T[lo] = Ti- {G'TlG'^Tr}. 



(2) 



(3) 



/(w) is the Bose-Einstein distribution function. The re- 
tarded Green's function C" is obtained from 

G'-H = {{u; + zO+)2 - K^' S£M - S'^H)"' , (4) 

where the retarded self-energy of lead a is 

i];M = y^"5^Hy"^, (5) 

and the lead surface Green's function g^[w] can be cal- 
culated from the generalized eigenvalue methodic, e.g., 
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kii and k^i are the diagonal and off-diagonal parts of 
the right lead dynamical matrix. Their sizes are deter- 
mined by the degrees of freedom M in each period of 
the lead. The matrix F^{s) translates to the right the 
displacement in the nth period to the (n -I- s)th period 
u]^{n+s) = F^{s)u]^{n). It is constructed from the eigen 
values and vectors of the generalized eigen value problem 



(u; + tO+)2/-A:f, 
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as F+is) = E+AXiEJ 



(7) 
. Here / is an M x AI iden- 
tity matrix. The diagonal matrix A^ consists of all the 
eigen values |A^| < 1, and E^ the corresponding eigen 
vectors E'l^ ~ (e^,e^, • ■ • ,£|//). Note that M' may be 
less than M, in which case the matrix inverse becomes 
pseudo-inverse. A similar left-translation matrix Fj^ (s) 
can be constructed from Ai and Ej^, which include all 
the eigen values |Ai| > 1 (excluding infinity) and the 
corresponding eigen vectors. 

While the NEGF method is systematic and suitable to 
take into account the nonlinear interaction, the transmis- 
sion coefficient from it is the sum of all the eigen modes 
from the leads. Thus it is difficult to analyse the con- 
tribution from each mode. The mode-matching method 
provides another way to calculate the transmission co- 
efficient in the ballistic limitiiiiSiii^. Single mode trans- 
mission and mode mixing effect can be studied by this 
method. The transmission matrix is given by 
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and 



RL 



{E^)-^g^RV^^G^V^^glkUFt{-l)-F^{-l))Et 

(9) 
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where ■;;+ ^ is the group velocity of the mth right prop- 
agating mode for the lead a. While the mode indices 
m and n are only for propagating modes, the matrices 
E'^ include all the propagating and evanescent modes. 
The total transmission coefficient as in the Landauer for- 
mula is T{u\ =. Tr{(t^^)tf«^}. Similar relations hold for 
waves incidcntcd from the right lead and transmitted to 
the left t^«. 

The NEGF and the mode-matching method are ex- 
actly equivalent in the ballistic case as show in Ref. [l^. 
All the matrices needed by Eqs. ([8H9]) can be obtained 
by solving the generalized eigen value problem Eq. ([7]). 
It is interesting to note that by doing a singular value de- 
composition on the transmission matrix t we can get the 
transmission eigenchannel information without any other 
efforts, for which different methods have been developed 
in the electronic transport literaturo^i. 
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FIG. 1: Ballistic thermal conductance of the benzene junction 
terminated by H atoms at different configurations. 3 is the 
full relaxed structure. 1 — 2 are in the compression state, and 
3 — 4 are in the stretching state. 



III. NUMERICAL RESULTS AND DISCUSSION 

We now present our numerical results. We begin with 
the simplest case where there is only one benzene ring at 
the center. We compare the thermal conductance of the 
molecular junction at different compression or stretching 
configurations by changing the distance between the two 
leads. The purpose of this study is twofold. Firstly, we 
want to study the effect of the compression or stretching 
on the thermal conductance. Secondly, our optimization 
process discussed in Sec. |TT] can not ensure that we have 
reached the lowest energy configuration while keeping the 
two SWCNTs to be coaxial. This is because we have fixed 
the position of the outer C and H atoms at both sides. If 
the thermal conductance is not sensitive to the distance 
between the two leads, our results make sense even if we 
do not find the lowest energy configuration. As we can 
see in Fig. [1] this is indeed the case. In Fig. [TJ from 1 to 
5 the distance between the two carbon atoms connecting 
to the amide groups is 10.07, 10.21, 10.38, 10.51, 10.64 A, 
respectively. The inset shows the atom configuration of 
the five cases. In the compression states (10.07 and 10.21 
A), the relative position between the benzene ring and 
the SWCNT leads changes compared to the full relaxed 
(10.38 A) and stretching (10.51 and 10.64 A) states (in- 
set of Fig. [1]) . Although the atom configurations change 
a lot, the ballistic thermal conductance of the molecu- 
lar junction only changes slightly in all the temperature 
range studied here. So we can conclude that the ballistic 
thermal conductance of the molecular junction is not sen- 
sitive to small compression or stretching of the molecule. 

In a related experiment^, the cutting edge carbon 
atoms are expected to be saturated by the COOH group, 
not the H group shown in Fig. [TJ In Ref. [23, the au- 
thors show that the electron transmission is largely af- 
fected by the terminating groups. It is also interesting 
to know how the terminating group affects the quantum 
thermal conductance. We still use the single benzene 
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FIG. 2: Ballistic thermal conductance of the benzene junc- 
tions terminated by H and COOH group. 



structure to study this problem. Figure [5] shows the ther- 
mal conductance of the two configurations. The full re- 
laxed structures are shown in the inset. The electrical 
conductance is mainly determined by the energy channel 
near the chemical potential, while the thermal conduc- 
tance is jointly contributed by many phonon modes. The 
terminating group only has a large influence on the high 
energy (short wavelength) optical phonon modes. The 
low energy (long wavelength) phonon modes are not sen- 
sitive to the local environment at the cutting edges. As a 
result, the thermal conductance only shows quantitative 
difference, which becoines larger at high temperatures. 

We now compare the conductance of the benzene rings 
and that of the alkane chains, with the SWCNTs termi- 
nated by H atoms. We expect that the COOH terminat- 
ing structures show similar behavior. To ensure the two 
molecules have comparable length, we include two ben- 
zene rings for the benzene structure and 8 CH2 groups 
for the alkane chain. The total transmission probability 
and the thermal conductance are shown in Figs. [3] and 31 
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FIG. 3: Phonon transmission probability as a function of en- 
ergy for the benzene and alkane chains. 
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FIG. 4: Ballistic thermal conductance as a function of tem- 
perature for the benzene and alkane chains. The inset shows 
the crossing point of the thermal conductance at about 38 K. 



respectively. Since the nonlinear interaction is not con- 
sidered here, phonon modes with different energies are 
independent of each other. The transmission probabil- 
ity is nonzero only in the overlapped energy range of 
the SWCNT leads and the center molecule. Strong cou- 
pling with the leads gives rise to wide broadening and 
strong shifting of the discrete molecule phonon energy. 
It is hard to make a one-to-one correspondence between 
the transmission peaks and the isolated molecule phonon 
eigen frequencies. This is especially true for low energy 
phonon modes, which have relatively large spatial ex- 
tent. The high energy modes are highly localized and 
show sharp peaks in the transmission spectrum. Due to 
the highly localized nature and low Bose-Einstein weight- 
ing, these modes can hardly transfer energy. There is a 
wide zero transmission gap around 0.07 eV for the alkane 
chain, which is a characteristic of the alkane chain phonon 
spectrum^. 

To gain further insight into the transmission spectrum, 
we also study the single mode transmission using the 
mode-matching method. Figure [5] shows the transmis- 
sion of several important branches in SWCNT phonon 



c 0.2 - 




0.03 0.04 0.01 

energy (eV) 



FIG. 5; Transmission probability of different phonon 
branches; (a) TA modes, (b) LA mode, (c) twisting mode, 
(d) lowest optical modes. 



spectrum: Two transverse acoustic (TA) phonon modes, 
one longitudinal acoustic (LA) mode, one twisting mode 
and the lowest two degenerate optical modes. These low 
frequency modes are expected to contribute much to the 
thermal conductance. The degenerate phonon modes of 
SWCNTs show different transmission in both cases. This 
is consistent with the fact that the junction structure de- 
stroys the degeneracy in transverse directions. At the en- 
ergy range below 0.01 eV, the main contribution comes 
from the TA modes of the SWCNT. The alkane chain 
shows larger transmission than the benzene chain in this 
energy range. Above 0.01 eV, the benzene chain shows 
larger transmission in most cases. The twisting mode 
shows the largest difference in two structures. While it 
can hardly be transmitted by the alkane chain, it has a 
large contribution in the benzene chain. 

The thermal conductance of the molecular junction de- 
pends not only on the transmission coefficient, but also 
on the phonon occupation number, which is reflected as 
the Bose-Einstein distribution in the Landauer formula. 
According the analysis of the transmission spectrum, we 
may expect that the alkane chain has larger thermal con- 
ductance at low temperatures, while at higher tempera- 
tures the order reverses. This is confirmed in Fig. |4l 
The crossing temperature is about 38 K. The room tem- 
perature thermal conductance is about 0.075 nW/K for 
alkane chains and 0.125 nW/K for benzene chains. In a 
recent experiment^, the thermal conductance of alkane 
chain is found to be smaller than our theoretical value. 
This difference comes from the effect of the leads. We 
are using SWCNTs here, while in the experiment it is 
the bulk gold connected with the alkane chain via the 
sufur atom. At least two factors from the leads may ac- 
count for this decrease. The first is the smaller phonon 
spectrum overlap between gold and alkane chain, and the 
second is the weaker coupling between them. 

Some comments are worthwhile. The nonlinear inter- 
action may change the transmission spectrum, and lead 



to a decrease of the thermal conductance. A perturbativc 
analysis of the single benzene structure shows that the 
room temperature thermal conductance drops about 30% 
of the ballistic value if we include the cubic force constant 
calculated from GaussionOS. Mean- filed approximation 
in the NEGF formalism only works well for simple struc- 
tures with relatively weak nonlinear interaction^. It fails 
to converge in the self-consistent iteration in present case. 
So it is still a challenge to find a good approximation for 
the nonlinear self-energies in NEGF method. 



NEGF or mode-matching method. Further information 
can be obtained from the transmission spectrum of each 
single mode. Using this method we show that the ballis- 
tic thermal conductance of the benzene ring amide linked 
with SWCNTs is not sensitive to the distance between 
the two SWCNTs. The benzene rings show larger ther- 
mal conductance than the alkane chains. This method 
is general, and can be easily applied to other material 
systems. 



IV. CONCLUSIONS 
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